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ABSTRACT 

We report the results of optical monitoring for a sample of 11 blazars including 10 
BL Lacs and 1 Flat Spectrum Radio Quasar (FSRQ). We have measured the multi- 
band optical flux and colour variations in these blazars on intra-day and short-term 
timescales of months and have limited data for 2 more blazars. These photometric ob- 
servations were made during 2009 to 2011, using six optical telescopes, four in Bulgaria, 
one in Greece and one in India. On short-term timescales we found significant flux 
variations in 9 of the sources and colour variations in 3 of them. Intra-day variability 
was detected on 6 nights for 2 sources out of the 18 nights and 4 sources for which we 
collected such data. These new optical observations of these blazars plus data from our 
previous published papers (for 3 more blazars) were used to analyze their spectral flux 
distributions in the optical frequency range. Our full sample for this purpose includes 
6 high-synchrotron- frequency-peaked BL Lacs (HSPs), 3 intermediate-synchrotron- 
frequency-peaked BL Lacs (ISPs) and 6 low-synchrotron-frequency-peaked BL Lacs 
(LSPs; including both BL Lacs and FSRQs). We also investigated the spectral slope 
variability and found that the average spectral slopes of LSPs show a good accor- 
dance with the Synchrotron Self-Compton (SSC) loss dominated model. Our analysis 
supports previous studies that found that the spectra of the HSPs and FSRQs have 
significant additional emission components. The spectra of all these HSPs and LSPs 
get flatter when they become brighter, while for FSRQs the opposite appears to hold. 
This supports the hypothesis that there is a significant thermal contribution to the 
optical spectrum for FSRQs. 

Key words: galaxies: active - BL Lacertae objects: general - galaxies: jets - galaxies: 
photometry 



1 INTRODUCTION 

Blazars are a very active class of active galactic nuclei, con- 
sisting of optically violently variable quasars, flat spectrum 
radio quasars (FSRQs) and BL Lacertae objects. In some 
cases, their emission lines are very weak or absent (equiva- 
lent width < 5 A ), a property which classically was used to 
define the BL Lacertae objects (Stocke et al. 1991; Marcha 
et al. 1996). Other blazars have broad emission lines similar 
to those of normal quasars. They have been observed at all 
wavelengths, from radio through very high energy (VHE) 
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7-rays. Blazars exhibit variability at all wavelengths on var- 
ious timescales. The high inferred isotropic luminosities and 
apparent superluminal motion revealed by long baseline ra- 
dio interferometry provide conclusive evidence that blazars 
are sources with relativistic jets oriented close to our line 
of sight and their emission is highly beamed and Doppler 
boosted in the forward direction. 

The overall emission, from radio to 7-rays shows the 
presence of two well-deflned broad components (von Mon- 
tigny et al. 1995; Fossati et al. 1998). The dominant emis- 
sion mechanisms in all classes of blazars is most likely syn- 
chrotron emission from radio to UV/soft X-ray frequencies 
and inverse compton scattering for hard X-ray and gamma 
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ray energies. The spectral energy distribution (SED) of LSPs 
(low synchrotron frequency peaked BL Lacs) and HSPs 
(high synchrotron frequency peaked BL Lacs) are system- 
atically different, with the observed peak of the emitted 
power typically at NIR/optical wavelengths for LSPs and 
at UV/soft X-ray wavelengths for HSPs (Giommi, Ansari 
& Micol 1995). One can define the objects as LSPs or 
HSPs if their ratio of X-ray flux in the 0.3-3.5 KeV band 
to radio flux density at 5 GHz is smaller than or larger 
than 10^^^'^, respectively (Padovani & Giommi 1996). How- 
ever, SED peaks can be located at intermediate frequen- 
cies as well, giving rise to the intermediate synchtrotron 
peaked blazar (ISP) classification (Sambruna, Maraschi & 
Urry 1996). Nieppola, Tornikoski & Valtaoja (2006) clas- 
sify over 300 BL Lacs and suggest that blazars with energy 
peak frequency, Vjj^ak ~ 10^^~^^ Hz are LSPs, those with 
Vpeak ~ lO^'^"^® Hz are ISPs and those with Vpeak ~ lO""^** 
Hz are HSPs. Recently, Abdo et al. (2010) extended the defi- 
nition to all types of non-thermal dominated AGNs depend- 
ing on the peak frequency of their synchrotron hump, Vpeak, 
and suggest the following classification: low-synchrotron- 
peaked (LSP) sources, with Upeak ^ 10^* Hz; intermediate- 
synchrotron peaked (ISP) sources, with 10^* <i^peak < 10^^ 
Hz; and high-synchrotron-peaked (HSP) sources, with i^peak 
^ 10^^ Hz; however, we have used the broader definition of 
ISPs put forward by Nieppola et al. (2006). 

Optical observations offer a wealth of information on the 
variability of blazars and play an important role in discrim- 
inating between LSPs and HSPs. Differences in the optical 
spectral properties of LSPs and HSPs constrain the theoret- 
ical models that try to explain the spectra and variability of 
these blazars. 

The study of SED properties is an excellent diagnostic 
tool for theoretical models and for the understanding of the 
physical radiation mechanisms, as these studies are crucial 
in analzing individual emission components. With respect 
to the optical spectrum of blazars, for HSPs, the maximum 
is located at higher frequencies and so the optical emission 
is in the ascending part of the synchrotron emission hump, 
while for LSPs the optical emission is expected to be in the 
descending part of this portion of the SED. Thus we expect 
for HSPs that the radio-through- X-ray spectral index orx is 
less than or equal to 0.75 or 0.80, while for LSPs, it is greater 
than 0.75 (Urry & Padovani 1995) or 0.80, whereas for ISPs, 
0.7^ urx ^0.8 (Sambruna, Maraschi & Urry 1996). 

Although the optical band is very narrow with respect 
to other spectral bands, it may yield a large amount of in- 
formation, as it can indicate the possible presence of other 
components in addition to the synchrotron continuum. For 
example, thermal emission from the accretion disk around 
the central engine is an important physical component, while 
the emission from the surrounding region of the nucleus 
and the host galaxy contribution are fundamentally contam- 
inants. The synchrotron inverse-Compton emission model 
predicts that the spectrum gets harder as the source turns 
brighter and several investigators indeed found that the op- 
tical spectrum became flatter when the flux increased, while 
the spectrum became steeper when the flux decreased (e.g.. 
Gear, Robson & Brown 1986; Maesano et al. 1997). How- 
ever, Ghosh et al. (2000) suggested that this might not be 
always correct. Trevese & Vagnetti (2002) analyzed the spec- 
tral slope variability of 42 PG quasars and concluded that 



the spectral variability must be intrinsic to the nuclear com- 
ponent. Vagnetti et al. (2003) showed that the spectral vari- 
ability, even restricted to the optical band, could be used 
to set limits on the relative contribution of the synchrotron 
component and the thermal component. 

Here we report the optical observations of 11 blazars on 
intra-day variability (IDV) and short-term variability (STV, 
from days through months) timescales during the 2009-2011 
observing seasons. We also searched for colour variations in 
these blazars on intra-day and short term timescales. In ad- 
dition, we analyzed the main features of the optical spectra 
for our sample of blazars. 

The paper is structured as follows. In Section 2, we 
present the observations and data reduction procedure, in- 
cluding our approaches to variability and time-scale detec- 
tion. Section 3 provides our results of IDV, STV and colour 
variability of our sample of blazars. Section 4 gives the gen- 
eral characteristics of the optical spectra. Section 5 provides 
the results of flux, colour and spectral variability of these 
sources and we have presented our conclusions in Section 6. 



2 OBSERVATIONS AND DATA REDUCTION 

Observations of these blazars were performed using six op- 
tical telescopes, four in Bulgaria, one in Greece and one in 
India. All of these telescopes are equipped with CCD de- 
tectors and Johnson UBV and Cousins RI filters. Details of 
five of these telescopes, detectors and other parameters re- 
lated to the observations are given in Table 1 of Gaur et al. 
(2012b) and a description of the remaining telescope is given 
in Table 1. A complete log of observations of these blazars 
from those six telescopes are given in Table 2. Although we 
have observed 13 blazars as part of this study, for our anal- 
yses of STV and colour variations, we used only 11 blazars, 
since the blazars J0211+1051 and lES 0806+524 were each 
only observed on 4 nights. 

We carried out optical photometric observations dur- 
ing the period 2009-2011. The raw photometric data were 
processed by standard methods. Image processing or pre- 
processing was done using standard routines in the Image 
Reduction and Analysis FacilitjQ (IRAF) and ES0-MIDAq3 
softwares. 

We processed the data using the Dominion Astrophys- 
ical Observatory Photometry (DAOPHOT II) software to 
perform the circular concentric aperture photometric tech- 
nique (Stetson 1987, 1992). For each night we carried out 
aperture photometry with four different aperture radii, i.e., 
IxFWHM, 2xFWHM, 3xFWHM and 4xFWHM. On com- 
paring the photometric results we found that aperture radii 
of 2xFWHM almost always provided the best S/N ratio, so 
we adopted that aperture for our final results. 

For these blazars, we observed three or more local stan- 
dard stars on the same fields. We employed two standard 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

^ ESO-MIDAS is the acronym for the European Southern Obser- 
vatory Munich Image Data Analysis System which is developed 
and maintained by European Southern Observatory 



Flux and spectral variations in blazars 3 



stars from each blazar field with magnitudes similar to those 
of the target and plotted their differential instrumental mag- 
nitudes for the IDV light curves. As the fluxes of the blazars 
and the standard stars were obtained simultaneously and so 
at the same air mass and with identical instrumental and 
weather conditions, the flux ratios are considered to be very 
reliable. Finally, to calibrate the photometry of the blazars, 
we used the one standard star that had a colour closer to 
that of the blazar. 

2.1 Vciriability Detection Criterion 

To search for and describe blazar variability wc have em- 
ployed two quantities commonly used in the literature. 

2.1.1 C-test 

The variability detection parameter, C was introduced by 
Romero et al. (1999), and is defined as the average of CI 
and C2 where 

Cl= -(BL-starA) ^ ^ ajBL - starB) 
a[starA — starB) ij{starA — starB) 

Here, {BL — star A) and {BL — starB) are the differential 
instrumental magnitudes of the blazar and standard star 
A and the blazar and standard star B, respectively, while 
a{BL — starA), a{BL — starB) and a{starA — starB) are the 
observational scatters of the differential instrumental mag- 
nitudes of the blazar and star A, the blazar and star B and 
starA and star B, respectively. If C > 2.576, the nominal 
confidence level of a variability detection is > 99%. However, 
this C-test is not a true statistic as it is not appropriately 
distributed and this criterion is usually too conservative (de 
Diego 2010). 

2.1.2 F-test 

We test our results for variability using the standard F-test, 

which is a properly distributed statistic (de Diego 2010). 
Given two sample variances such as Sq for the blazar in- 
strumental light curve measurements and s* for that of the 
standard star, then 



The number of degrees of freedom for each sample, pq 
and will be the same and equal to the number of mea- 
surements, N, minus 1 (i^ = N — 1). The F value is then 
compared with the -Fig'^^ critical value, where a is the sig- 
nificance level set for the test. The smaller the a value, the 
more improbable that the result is produced by chance. If 
F is larger than the critical value, the null hypothesis (no 
variability) is discarded. We have performed f-tests at two 
significance levels (0.1% and 1%) which correspond to 3cr 
and 2.6cr detections, respectively. 

2.2 Time-Scale of Variability Detection Methods 

2.2.1 Structure Function 

The structure function (SF) is a useful tool to make prelim- 
inary estimates of any periodicities and timescales in time- 



series data and is able to discern the range of the charac- 
teristic timescales that contribute to the fluctuations. For 
details about how we used the SF see Gaur et al. (2010). 
Emmanoulopoulos et al. (2010) have discussed the weak- 
nesses of the SF method, so wc have cross checked our SF 
results by the discrete correlation function (DCF) method. 

2.2.2 Discrete Correlation Function 

The DCF is analogous to the classical correlation function 
which requires evenly sampled data except that it can work 
with unevenly sampled data. It was introduced by Edelson 
& Krolik (1988). Hufnagcl & Brcgman (1992) gencrahzcd 
the method to include a better error estimate. For details 
about the DCF see Tonnikoski et al. (1994), Hovatta et al. 
(2007) and references therein. 



3 RESULTS 

3.1 Intra-day, Short-term flux and Colour 
veiriability of individual bleizeirs 

3.1.1 3C 66A: 

3C 66A was classified as a BL Lac object by Maccagni et al. 
(1987), based on its significant optical and X-ray variability 
with its redshift at z = 0.444 (Lanzetta et al. 1993) that 
value is actually quite uncertain (e.g., Bramcl et al. 2005). 
The synchrotron peak of this source is located between 10^^ 
and 10^® (Perri et al. 2003) and therefore 3C 66A can be 
classified as an ISR It shows highly variable non-thermal 
continuum emission ranging from radio up to X-ray and to 
7-ray frequencies on different timescales and exhibits strong 
polarisation from radio to optical wavelengths (e.g. Maraschi 
et al. 1994; Urry & Padovani 1995; Hartman et al. 1999; 
Gupta et al. 2008; Villata ct al. 2009; Bauer ct al. 2009). 
Long term optical variability and spectral variability have 
been studied before (Fan & Lin 2000; Vagnetti et al. 2003; 
Hu et al. 2006). Short-term variability and color variations 
have been investigated by many authors (e.g. Takalo et al. 
1996; Gu et al. 2006; Rani ct al. 2010a). 

We observed 3C 66A for STV in 16 nights between 13 
November 2009 to 6 November 2010. The source showed 
large flux variations, but no signiflcant colour variations dur- 
ing our observing run (Fig 1). The source showed a bright- 
ening trend by ~ 0.5 mag in 38 days in the R band and 
again was as bright as 13.4 mag even after another 319 days 
(around JD 2455505). This flux is comparable to the bright- 
est R magnitude (~13.4) reported in the source by Bottcher 
et al. (2009). So we conclude that 3C 66A was in a high 
state during this period. 

3.1.2 AO 0235+164: 

AO 0235-1-164 has a redshift oi z = 0.94 based on the de- 
tection of emission lines (Nilsson et al. 1996) and was clas- 
sified as a BL Lac object by Spinrad & Smith (1975). The 
synchrotron peak of this source is located at 10^^'^ Hz and 
therefore it is classified as an LSP (Nieppola et al. 2006). 
This blazar has exhibited extreme variability (by over an 
order of magnitude on time scales < 1 year) across all spec- 
tral ranges, including X-ray (Raiteri et al. 2009) and 7-ray 
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(i.e, Abdo et al. 2010). It has been observed from the radio 
to X-ray bands on timescales ranging from less than an hour 
to many years (e.g. Heidt & Wagner 1996; Fan &: Lin 1999; 
Romero et al. 2000; Webb et al. 2000; Raiteri et al. 2001; 
Padovani et al. 2004). Raitori ot al. (2001) predicted that 
the blazar should show a possible correlated periodic radio 
and optical outburst on the timescale of 5.7±0.5 years that 
was expected in 2004 February-March but that predicted 
outburst could not be detected (Raiteri et al. 2005, 2006a). 
Raiteri ct al. (2006b) then reanalyzed the long-term optical 
light curves and suggested that the optical outbursts may 
have a longer timescale of ~8 years. Gupta et al. (2008) re- 
ported an outburst in .January 2007, with Kmag — 14.98 the 
peak seen during their observing run. 

Wc have observed AO 0235+164 on 8 nights between 13 
November 2009 to 6 November 2010. and it showed signifi- 
cant flux as well as colour variations (Fig 1). This source's 
magnitude varied between 17.9 to 18.3 mag in R band, and 
then was at ~ 18.4 mag after ~347 days (at JD 2455508) 
which is ~4 mag fainter than the brightest reported value for 
this source. We conclude that AO 0235-1-164 was probably 
in a low state in the period we observed it. 

3.1.3 S5 0716+714: 

The blazar S5 0716-1-714 at redshift z = 0.31 ± 0.08 is one 
of the most active sources in optical bands. Its optical UV 
continuum is so featureless (Biermann et al. 1981; Stickel 
et al. 1993) that a redshift estimate has been possible only 
by resolving and using the host galaxy as a standard can- 
dle (Nilsson ot al. 2008). This BL Lac is of the ISP type, 
displaying a first broad peak in the optical/UV bands and 
showing another broad peak near 1 GeV (Massaro et al. 
2008; Ferrero et al. 2006). It is highly variable on timescales 
from hours to months across all observed wavelengths (e.g., 
Wagner et al. 1990; Heidt & Wagner 1996; Villata et al. 
2000; Raiteri et al. 2003; Montagni et al. 2006; Foschini et 
al. 2006; Ostorero et al. 2006; Gupta et al. 2008, and refer- 
ences therein). This source is one of the brightest BL Lacs 
in optical bands and has an IDV duty cycle of nearly 1 (e.g., 
Wagner et al. 1996; Stalin et al. 2006; Montagni et al. 2006). 
Gupta, Srivastava & Wiita (2009) analysed the excellent in- 
traday optical LCs of the source obtained by Montagni et 
al. (2006) and reported some evidence for nearly periodic 
oscillations ranging between 25 and 73 minutes on several 
different nights. More recently, Rani et al. (2010b) reported 
very good evidence for a quasi-periodic oscillation of ~15 
minutes in a new optical light curve of S5 0716-1-714. 

We observed S5 0716-1-714 on 19 nights between 14 
November 2009 to 25 March 2011 for short-term variability 
and in three nights in R pass-band (11, 13 and 15 December 
2010) for intra-day variability. This source exhibited signifi- 
cant STV and colour variations (Fig. 1). Also, the IDV light 
curves of the blazar and the differential instrumental magni- 
tude of standard stars with arbitrary offset are displayed in 
Fig. 2. The source showed genuine IDV on all three nights 
and the visual impressions of the LCs were confirmed by C- 
and F-tests for those nights (Table 3). We performed SF 
and DCF analyses to find if any time scale of variability 
was present in these nights (shown in Fig. 2) but no sig- 
nificant time scale was found. There have been four major 
optical outbursts reported from S5 0716-1-714 at the begin- 



ning of 1995, in late 1997, in the fall of 2001, and in 2004 
March (Raiteri et al. 2003; Foschini et al. 2006). These four 
outbursts give a possible period of long-term variability of 
~3.0±0.3 years. During our observations, the source first 
brightened by ^^^1.15 mag in 32 days in the R band and then 
it fades by ^1 mag over 27 days (covering the period JD 
2555149 to 2455231). Later on, within 94 days this source 
faded by ~0.98 mag (ending at JD 2455596) and then again 
it brightened and reached the highest flux we saw at 12.54 
mag at JD 2455625. This is essentially as bright as ever 
reported in the earlier outbursts in the R passband (12.55 
mag in the fall of 2001; Raiteri et al. 2003) and 12.55 mag in 
January 2007 (Gupta ct al. 2008). Wc can conclude that the 
source was in an outburst state during our observing run. 

3.1.4 OJ 287: 

OJ 287 is a LSP at redshift z = 0.306 (Sitko & Junkkarinen 
1985) and has received a lot of attention as it has shown 
massive double-peaked outbursts approximately every 12 
years since the 1890s. Sillanpaa et al. (1988) first noticed 
this exceptional behaviour and suggested that these out- 
bursts might be caused by a close binary black hole system 
in which the secondary black hole induces tidal disturbances 
in the accretion disc of the primary black hole. Several other 
authors suggested models for OJ 287 (Katz 1997; Villata ct 
al. 1998; Valtaoja et al. 2000), all of them based on the as- 
sumption that OJ 287 hosts a close binary black hole. It has 
been observed extensively in optical band (e.g., Carini et al. 
1992; Sillanpaa et al. 1996a,b; Abraham et al. 2000; Gupta 
et al. 2008; Fan & Lin et al. 2009). The observational prop- 
erties of OJ 287 from radio to X-ray energy bands have been 
reviewed by Takalo et al. (1994). Fan et al. (2009) reported 
large variations in the source with AV = 1.96 mag, AR = 
2.36 mag and Al = 1.95 mag during their observations span- 
ning 2002 to 2007. This quasar shows quasi-periodic optical 
outbursts at ~12 yr intervals, with two outburst peaks per 
interval. Valtonen et al. (2008) discussed a model in which 
a secondary body (a black hole) pierces the accretion disk 
of the primary black hole and produces two impact fiashes 
per period. This model predicted the next outbursts to occur 
between 2007 October and 2009 December and they were ac- 
tually seen (Valtonen & Ciprini 2011 and refernces therein) 

We observed OJ 287 on 16 nights between 13 Novem- 
ber 2009 and 2 November 2010 for STV and on five nights 
between 2009 November 22 and 2010 January 20 for IDV. 
OJ 287 showed significant fiux variations in all the observed 
pass-bands but no significant colour variations were found on 
short time scales (Fig. 1). Also, the IDV light curves of the 
source and the differential instrumental magnitude of two 
standard stars with arbitrary offset are displayed in Fig. 3. 
We performed C and F-tests and found that the source did 
not show genuine IDV on any of the five nights. Real STV 
was certainly seen, as the source brightened by ^ 0.6 mag 
in 38 days and peaked upto 13.56 mag in R band. Then, it 
faded to 15.01 mag in 317 days during our observing run. 
Pursimo et al. (2000) carried out multi-band optical mon- 
itoring of the source for about five years (1993-1998) and 
also plotted its V passband historical light curve for about 
century long observations. The 1993-1998 light curve of the 
source shows the brightest (~13.5) and faintest (~16.5 mag) 
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states of the source in the R pass-band. In our present ob- 
servations, we observed the source at R ~ 13.56 mag on JD 
2455187 in the R pass-band which is near to its brightest 
state. So, near December 2009, the source was in outburst 
state and then faded by ~ 1.54 mag by November 2010. 

3.1.5 Mrk 421: 

Mrk 421 has z = 0.031 and is classified as an HSP because 
the energy of its synchrotron peak is higher than 0.1 KeV. 
It is the brightest TeV 7-ray emitting blazar in the north- 
ern hemisphere. It was first noticed to be an object with a 
blue excess which turned out to be an elliptical galaxy with 
a bright point like nucleus (Ulrich et al. 1975). It was the 
first known extragalactic TeV 7-ray emitter (Punch et al. 
1992). It has been extensively observed at all wavelengths 
from radio to 7-rays. The source is characterized by strong 
variability in the optical region (i.e., Miller 1975; Liu et al. 
1997) including optical variations of at least 4.6 mag (Stein 
et al. 1976). Mrk 421 has been a target of several simulta- 
neous multi-wavolcngth campaigns (Takahashi ct al. 2000; 
Rebillot et al. 2006; Fossati et al. 2008; Lichti et al. 2008). 

We observed Mrk 421 for STV on 17 nights between 15 
November 2009 and 2 May 2011 in B, V, R and I bands. 
The flux from the nucleus of Mrk 421 is contaminated by 
the omission of the host galaxy so wc used the measure- 
ments of Nilsson et al. (2007) to estimate the host galaxy 
emission in the R-band. This flux is used to obtain the cor- 
responding contribution for B, V and I bands (Fukugita et 
al. 1995) and we corrected for the Galactic extinction using 
the extinction map of Schlegel, Finkbeiner & Davis (1998). 
The source showed genuine flux variations in the V, R and 
I bands, but no significant colour variations during our ob- 
serving run (Fig. 1). During our observing run, the source 
showed an average brightness of ~13.4 mag and reached a 
maximum brightness of ^12.46 mag in the V band at around 
JD 2455684 which is quite close to the maximum observed 
brightness of about V'^12.0 in November 1996 during an 
outburst (Tosti ct al. 1998a). Thus we can say that this 
source was in an outburst state during our observing run. 

3.1.6 lES 1218+304: 

lES 1218-1-304 has a redshift of 0.182 and was first detected 
to emit VHE 7-rays by MAGIC in 2005 (Albert et al. 2006). 

It is classified as an HSP as its synchrotron peak is located 
at 10" " Hz (Nicppola ct al. 2006). In 2009 VERITAS re- 
ported fgist variability from the source, the peak flux reach- 
ing ~26 per cent of the Crab Nebula Flux (Acciari et al. 
2010). 

We observed the source on 6 nights between 2 March 
2011 and 7 June 2011 (Fig. 1). We performed C- and F- 
tests reported in Table 4 on these nights but neither STV 
nor colour variations were found for this source during our 
observations. 

3.1.7 ON 231: 

ON 231 (also known as W Comae) was discovered in 1971 
because of its odd properties, particularly a strong opti- 
cal variability and an apparently lineless continuum (Bi- 



raud 1971; Browne 1971). The detection of weak emission 
lines in its spectrum made it possible to determine that its 
z = 0.102 (Weistrop ct al. 1985). The historial B passband 
LC of the source (during 1935-1997) was plotted by Tosti et 
al. (1998b). It has shown optical flux variations on diverse 
timescales ranging from a few hours to several years (Xie 
et al. 1992; Smith & Nair 1995). After the optical outburst, 
ON 231 showed a slow decline in its mean luminosity (Tosti 
et al. 2002). This source was also detected at high energies 
by Acciari ct al. (2008). 

We observed ON 231 for STV on 17 nights between 15 
November 2009 and 6 February 2011 in B, V, R and I bands. 
The light curves of the calibrated blazar and differential in- 
strumental magnitude of two standard stars with arbitrary 
offset are displayed in Figure 4. During our observations ON 
231 showed significant flux, but no signiflcant colour, vari- 
ations. Tosti et al. (1998b) carried out multi-band optical 
monitoring of the source for three years (1994-1997) in its 
great outburst state and gave their brightest (~ 13.5 mag) 
and faintest (~ 15.0 mag) values in the R passband. In the 
present observations, the source brightened by ^0.3 mag in 
the R pass-band over ~11 days and peaked at 14.41 mag 
which is ~0.6 mag fainter than the faintest state observed 
by the Tosti et al. (1998b). Then, the source faded to -15.08 
mag at JD 2455600. It appears likely we have observed ON 
231 in a low state. 



3.1.8 3C 279: 

The FSRQ 3C 279 shows flux variabilities at all frequencies. 
It is a LSP having synchrotron peak frequency at 10^^'* 
Hz (Planck CoUaboration et al. 2011). Webb et al. (1990) 
reported rapid fluctuations of ~2 mag within 24 hours at 
visible wavelengths and a rapid variation of 1.17 mag within 
40 min was seen in the V band (Xie et al. 1999). Gupta et al 
(2008) reported a 1.5 mag variation in the R band over 42 
days and Rani et al. (2010a) reported a rapid decay in the 
brightness of the source in 5 days by AR=1.10 mag. This 
source has been studied through multiwavcleiigth campaigns 
(Hartman et al. 1996; Wehrle et al. 1998; Bottcher et al. 
2007). 

Wc observed 3C 279 on 11 nights between 27 May 2011 
and 24 June 2011. We performed C- and F-tests on these 
nights and found genuine STV in V, R and I bands shown 
in Fig 4. The source did not show genuine colour variations 
during the observations. In the R band, 3C 279 first faded 
by ~0.47 mag in 16 days from 14.36 mag to 14.83 mag and 
then brightened again to 14.62 mag in 15 days. The source 
was reported to in outburst in 2007 January by Gupta et 
al. (2008) reaching a brightness of R ~ 12.6 mag and Rani 
et al. (2010a) has reported 17.1 mag as the faintest level in 
the source. So wc can say that 3C 279 was in either a pre- 
or post-outburst state during our observations. 

3.1.9 lES 1426+428: 

The blazar lES 1426-1-428 was discovered in the medium 
X-ray band (2-6 KcV) with the Large Area Sky Survey Ex- 
periment (LASS) on the first of the High Energy Astron- 
omy Observatories (Wood et al. 1984). The host galaxy of 
lES 1426-1-426 was resolved by Urry et al. (2000) and found 
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to have mii=16.14, leading to a photometric redshift of 
z=0.132±0.030 It was classified as a BL Lac object because 
of its lack of prominent optical emission lines by Remillard 
et al. (1989). Based on observations by BeppoSAX upto 
100 KeV, lES 1426+428 is an example of an extreme-high- 
energy peaked BL Lac object (HSP) with the peak of the 
synchrotron emission lying above 100 KeV (Costamante et 
al. 2001) during that observing campaign. Hence it was a 
prime candidate for TeV gamma-ray emission, which was 
confirmed by the Whipple Collaborations (Horan et al. 2002; 
Petry et al. 2000). Optical polarization and photometry of 
the source was done by Jannuzi et al. (1993, 1994). 

We observed lES 1426-1-428 for short-term variability 
on 14 nights between 16 March 2010 to 8 February 2011 
out of which we observed the source for intra-day variability 
on four nights (15, 17 March 2009 and 18, 21 June 2010). 
We found genuine STV and colour variation in this source. 
The light curve of the source for STV and its colour vari- 
ations are shown in Fig. 4. IDV light curves of the blazar 
and the differential instrumental magnitudes are displayed 
in Fig. 2. For two nights (18 and 21 June 2010), we have 
observed the source quasi-simultaneously in the B and R 
bands. The source showed genuine IDV during three nights 
which are confirmed by C— and F-tests for these nights. We 
performed SF and DCF analyses to search for any timescales 
of variability but no significant timescale was found. During 
our observing run, the source showed an average magnitude 
of ~15.9 in the R band. In last decade, the source has only 
varied between ~15. 7-16.1 mag[f|. So, the source in our ob- 
servations also shows a similar state of brightness. 

3.1.10 lES 1553+113: 

IBS 1553+113 was discovered in The Palomar-Green Survey 
of ultraviolet-excess objects as a 15.5 magnitude blue stellar 
object (Green, Schmidt & Liebert 1986) and BL Lac clas- 
sification was suggested by the featureless spectrum (Miller 
& Green 1983) with an optical R magnitude varying from 
~13 to ~15.5 (Miller et al. 1988). Its SED indicates that 
it is an HSP (Falomo & Treves 1990; Donoto et al. 2005). 
The logarithmic ratio of its 5 GHz radio flux, Fsghz, to 
its 2 keV X-ray flux, F2KeV, has been found to range from 
log{F2Kev/F5GHz)^ -4.99 to -3.88 (Osterman et al. 2006; 
Rector et al. 2003). lES 1553+113 has been detected from 
radio through hard X-rays and also in the very high en- 
ergy (VHE; E ^ 100 GeV) bands up to energies above 1 
TeV (Aharonian et al. 2006; Albert et al. 2007). It is in the 
Fermi LAT bright active galactic nucleus (AGN) source list 
(LBAS; Abdo et al. 2009). 

It is a bright optical source with V-band magnitude 
of V ~ 14 (Falomo & Treves 1990; Osterman et al. 2006). 
Observations taken between 1986 and 1991 found its optical 
spectral index, ao, to remain constant (ao ~ — 1) and its 
magnitude to vary by AV=1.4 (Falomo et al. 1994). 

We observed lES 1553+113 for STV on 15 nights be- 
tween 18 March 2010 and 6 August 2010 and in 6 nights for 
IDV (between 2010 May 7 and 2010 June 21). In four out of 
these six nights, we observed the source in B and R bands 
quasi-simultaneously and also searched for colour variations 

^ http://users.utu.fi/kani/lm/ 



in (B-R) band on intra-day time scales (Fig. 3). We per- 
formed C and _F-tests in all nights but genuine IDV was not 
found during any night. Also, colour variations in (B-R) were 
not found on these four nights. Genuine STV was found but 
the source exhibited no significant colour variations (Fig. 4) . 
During our observing run, the source exhibited the average 
magnitude of ~13.3 mag in R band. This value is even higher 
than the brightest earlier measurement (13.5 mag) observed 
by Osterman et al. (2006) during the flaring state of source. 
So we clearly have observed this source in a high state. 



3.1.11 BL Lac: 

BL Lacertae (2200+420), lies in an elliptical galaxy at a 
redshift of -0.07 (Miller, French & Hawley 1978). It is clas- 
sified as a LSP (Fossati et al. 1998). BL Lacertae has been 
studied in detail during various intensive multiwavelength 
campaigns (Madejski et al. 1999; Ravasio et al. 2003; Vil- 
lata et al. 2002; Bottcher et al. 2003 and references therein). 
Fan et al. (1998) presented historical optical data (in UB- 
VRI) over ~100 years covering the period 1896-1996 and 
found large variations in all the bands. 

We observed this prototypical source on 6 nights be- 
tween 15 July 2009 and 17 June 2010. BL Lac exhibited 
significant flux variations but no colour variations during 
our observing period. In the R band, the source was in a 
low state at 14.83 mag on JD 2455027 which is quite close 
to the faintest value (14.91 mag) mentioned in Fan & Lin 
(2000). During our observations on 5 days (12-17 Jun 2010) 
source showed the average magnitude of ~ 13.40 mag in R 
band which is only ~1.5 mag brighter than the faintest mag 
reported for the source by Fan & Lin (2000). So it is fair to 
say that the source was in a low state during our observing 
run. 



3.2 Correlated variations between colour and 
magnitude 

We have seen the corresponding variations in the (B-V), (V- 
R), (R-I) and (B-I) colour indices of these blazars with re- 
spect to variations in brightness. Colour-magnitude plots of 
the individual sources are displayed in Fig. 5. The individual 
panels show the (B-V), (V-R), (R-I) and (B-I) colours plot- 
ted (in sequence from bottom to top with arbitrary offsets) 
with respect to V magnitude. The straight lines shown are 
the best linear flt for each of the colour indices, Y , against 
magnitude, V, for each of the sources: Y — mV + c. Fitted 
values for the slopes of the curves, m, and the constants, 
c, are listed in Table 5 and 6. We have also quoted the lin- 
ear Pearson correlation coefficients, r, and the corresponding 
null hypothesis probability values, p. 

A positive correlation here is defined as a positive slope 
between the colour index and apparent magnitude of the 
blazar. This means that the source tends to be bluer when 
it brightens or redder when it dims. An opposite correlation 
(with negative slope) implies the source follows a redder 
when brighter behaviour. 

We found significant positive correlations {p ^0.05) be- 
tween the V magnitude and some of the color indices for the 
blazars: S5 0716+714 (B-V); OJ 287 (B-V); Mrk 421 (B-V), 
(R-I), (B-I), ON 231 (V-R), (R-I), (B-I) and lES 1426+428 
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(V-R), (B-I). AO 0235+164 (V-R) showed a weak positive 
correlation. No significant negative correlations are found 
for any source. 3C 66A, lES 1218+304, 30 279 and BL Lac 
lacked any significant correlations greater than 95 per cent 
confidence. For lES 1553+113, we found variations of (B-R) 
colour with respect to V magnitude but the linear Pearson 
correlation coefficient, r = —0.146 gives a p = 0.441; hence, 
no significant correlation is found for this source. 



4 THE OPTICAL SPECTRAL SLOPE 
DISTRIBUTION 

We have combined the optical data in the B, V, R and I 
bands from our present observations with optical observa- 
tions of three more sources: 2 HSPs (lES 1959+650 and 
lES 2344+514) and one LSP (30 454.3) from our papers 
Gaur et al. (2012a) and Gaur et al. (2012b), respectively. 
We can not use lES 1553+113 in this aspect of our study 
as the standard stars in its field arc available only in the B 
and R bands. We regard these quasi-simulataneous multi- 
band observations taken within 30 minutes as "simultane- 
ous" observations (Fiorucci et al. 2004), but we caution that 
occasionally a significant variation will have occurred be- 
tween the measurements at different bands. All the obser- 
vational data in the B, V, R and I bands were corrected 
for foreground Galactic interstellar reddening, which is de- 
duced from the maps of dust infrared emission reported by 
Schlegel, Finkbeiner & Davis (1998). 

The radiation of these blazars can be adequately de- 
scribed by a single power law Fi, = Av~" in the optical 
band, where a is the spectral slope which is obtained by us- 
ing the linear least square fitting between log and log v. In 
no case do we measure significant curvature to the spectra. 
We have accepted only the data that allow us to calculate 
the regression based on at least three quasi-simultaneous 
photometric bands, and whose square of Pearson's linear 
correlation coefficient is greater than 0.9. Also, the standard 
deviation of the fitted slope must be less than 0.25 for us 
to consider the fit valid. These selection criteria assure the 
accuracy of the slopes by cutting out all the sets of three or 
four data points that contain at least one wrong value (Hu 
et al. 2006). These results arc listed in Table 7. Our sam- 
ple consists of 6 LSPs, 3 ISPs and 6 HSPs. The sequence 
of columns in Table 7 are: the source name; each photomet- 
ric band average flux density (in mjy) and the number of 
useful data points; the average optical spectral index a, the 
standard deviation of a and the number of data points; the 
variation amplitude (max-min) of the optical spectral slope 
Ma ; the slope of the linear regression between a and the R 
magnitude, 6, and its standard deviation; the correlation co- 
efficient between a and R magnitude; and the corresponding 
p values. 



5 DISCUSSION 
5.1 Flux Variations 

5.1.1 Intra-Day and Short-term Variability Time-scales 

We carried out photometric observations of four blazars S5 
0716+714, OJ 287, lES 1426+428 and lES 1553+113 for 



IDV in the R band. For some of these nights, we have car- 
ried out quasi-simultaneous observations in B and R pass- 
bands for lES 1426+428 and lES 1553+113. IDV was not 
detected for OJ 287 and lES 1553+113 during any night as 
seen in Fig. 3. For two other sources, S5 0716+714 and lES 
1426+428, we found genuine IDV in 2 (out of 3) and in 4 
(out of 5) nights, respectively (shown in Fig. 2). We searched 
for time-scales of variability using SF and DCF techniques, 
but no significant variability of time-scale was found for any 
of these blazars. 

We performed photometric observations of nine blazars 
in the BVRI pass-bands during our observing run cover- 
ing the period 2009 2011. Comparisons of our observations 
with earlier measurements of the same sources indicated that 
the blazars 30 66A, S5 0716+714, OJ 287, Mrk 421 and 
lES 1553+113 were probably in high states; blazars AO 
0235+164, ON 231 and BL Lac were apparently in faint 
states and 30 279 was in a pre- or post-outburst state. 
The substantial flux variations for optical IDV and STV 
in blazars can be reasonably explained by models involving 
relativistic shocks propagating outwards (e.g. Marscher & 
Gear 1985; Wagner & Witzel 1995; Marscher et al. 1996), 
although most of the IDV in the radio band comes from 
interstellar scattering. The larger flares are expected to be 
produced by the emergence and motion of a new shock trig- 
gered by some strong variations in a physical quantity such 
as velocity, electron density or magnetic field moving into 
and through the relativistic jet. IDV reported in the blazars 
is generally attributed to the shock moving down the in- 
homogeneous medium in the jet. In AO 0235+164, gravi- 
tational microlensing (e.g. Gopal-Krishna & Subramanian 
1991) may play a role as it is known to have two galaxies 
(at z = 0.524 and 0.851) along our line of sight to it. 

5.2 Colour Variations 

We have reported searches for variations in colour with 
time for time periods corresponding to STV in a sample of 
blazars. Out of eleven blazars we have newly observed, only 
three (AO 0235+164, S5 0716+714 and lES 1426+421) have 
shown significant colour variations on short timcscales. Phe- 
nomena that could lead to colour variations with time were 
summarized by Hawkins (2002) who noted they could arise 
from: colour changes in accretion discs; different contribu- 
tions from the underlying host galaxy arising from changes 
in seeing; and colour changes from microlensing. 

Colour variations in these sources are most unlikely to 
arise from the underlying host galaxy, as in the case of BL 
Lacs, Doppler boosed jet emission almost invariably swamps 
the light from the host galaxy. Also, BL Lacs seem to have 
much weaker disc emission, so colour changes arising in the 
accretion disc are also unlikely. AO 0235+164 showed sig- 
nificant colour variations and is known to have two fore- 
ground galaxies at z — 0.524 and z — 0.851 (Nilsson et al. 
1996) so colour variations seen in this source could arise 
from microlensing of different regions within a relativistic 
jet (Gopal-Krishna & Subramanian 1991). In the case of BL 
Lacs, any accretion disc radiation is always overwhelmed by 
that from the jet, and therefore jet models that can produce 
different fluctuations in different colours are most likely to 
be able to explain the colour variations we have detected. 
In the usual models involving shocks propagating down the 



8 Gaur et al. 



jets (e.g. Marscher & Gear 1985; Marscher et al. 2008), ra- 
diation at different frequencies is produced at different dis- 
tances behind the shocks. Higher frequency photons from 
the synchrotron component emerge sooner and closer to the 
shock front. Hence colour variations in blazars could arise 
from differences in the times and rates at which radiation 
observed at different visible colours rise and fall. 

5.3 Relation between colour index and flux 
Vciriation 

We now investigate the relationship between spectral 
changes and flux variations on short-term as well as intra- 
day time scales. On short time scales, out of the eleven 
blazars for which new data is presented here, five (S5 
0716-F714, OJ 287, Mrk 421, ON 231 and lES 1426-^428) 
have shown significant positive correlation between colour 
and magnitude; one (AO 0235+164) has shown a very weak 
positive correlation, while five (3C 66A, lES 1218-1-304, 3C 
279, BL Lac & lES 1553-1-113) lacked any correlation. 

No significant negative correlation was found for any BL 
Lac which reflects the common trend for this class that they 
become bluer when they brighten (GhiscUini et al. 1997; Fan 
et al. 1998: Massaro et al. 1998; Fan & Lin 1999; Ghosh et 
al. 2000; Raiteri et al. 2001; Villata et al. 2002; Gu et al. 
2006; Rani et al. 2010a). It has been found that the ampli- 
tude of the variations is systematically larger at higher fre- 
quencies, which suggests that the spectrum becomes steeper 
when the source brightness decreases, and flatter when it 
increases (Ghisellini et al. 1997; Massaro et al. 1998). The 
variations of the colour indices of our BL Lac objects follow 
this trend. It is possible to explain this colour variation with 
a one-component synchrotron model, i.e., the more intense 
the energy release, the higher the typical particle's energy 
and the higher the corresponding frequency (Fiorucci et al. 
2004). It could also be explained if the luminosity increase 
was due to the injection of fresh electrons with an energy dis- 
tribution harder than that of the previous, partially cooled 
ones (e.g. Kirk et al. 1998; Mastichiadis & Kirk 2002). How- 
ever, some evidence that the amplitudes of variations are 
not systematically larger at high frequencies has been found 
on several occasions (Brown et al. 1989; Massaro et al. 1998; 
Ghosh et al. 2000; Ramirez et al. 2004). Ghosh et al. (2000) 
suggested that it may not be correct to generalize that the 
amplitude of the variation in blazars is systematically larger 
at higher frequency as they found a reddening as the bright- 
ness increases in their optical observations of the BL Lac 
objects PKS 0735-hl78, 3C 66A and AO 0235-M64. 

The trend of steeper spectra when brighter seems to 
generally hold for the variations we saw in IDV. On intra- 
day timescales, lES 1426+428 and lES 1553+113 (during 
one night) showed significant negative correlations between 
colour index and magnitude. So we can say that the trend 
of flatter spectra when brighter does not always hold for 
variations on time scales of days. While amplitudes of vari- 
ations are usually larger in the B band than in the R band, 
sometimes those amplitudes of variation were comparable. 

5.4 Spectral Variations 

In figure 6 we have plotted the spectral indices a in the opti- 
cal band against synchrotron peak frequency (vpeak) for all 



classes of blazars. The values of a have been taken from our 
new measurments of these sources and have been supple- 
mented by data for some more sources from Fiorucci et al. 
(2004) and Hu et al. (2006). Synchrotron peak frequencies 
of these sources has been taken from Nieppola et al. (2006) 
and Planck Collaboration et al. (2011). Spectral variations 
of these classes of blazars are discussed below. 

5.4.1 LSPs 
BL Lacs 

For LSPs, wc expect optical spectral indices aave >1 
due to their being on the descending part of the first bump 
of the SED. In this paxt of the spectrum we expect the 
spectral index to have broken so that a « athin + 1/2, 
where athin — 0.5 -0.7. Wc have 11 LSPs in this sample 
and indeed all save 1 have aave > 1- From figure 6, it is 
clear that the distribution of spectral indices in the optical 
clustered around 1.5. The average of the spectral slopes of 
these LSPs is ~1.5 excluding the two anomalous sources 
PKS 2032+107 and AO 0235+164. Chiang & Bottcher 
(2002) demonstrate that for a broad range of particle 
injection distributions, SSC-loss-dominated synchrotron 
emission exhibits exactly this kind of spectral slope. Hence, 
SSC models fit the optical emission for LSPs very well. 

FSRQs 

Since FSRQs are a subclass of LSPs, we expect the op- 
tical spectral index also should be greater than 1. However, 
FSRQs show strong emission lines and a thermal contribu- 
tion that may be comparable to the synchrotron emission in 
the optical spectral region. The optical emission of FSRQs 
is often strongly contaminated by thermal emission from the 
accretion disk and it is shown by the presence of the "big 
blue bump" in the optical/UV region. This contribution was 
clearly shown in the SED of 3C 273 by Tiirler et al. (1999) 
and the mean optical spectral slope was calculated to be 0.6 
by Fiorucci et al. (2004). However, for other two sources, 
3C 454.3 and 3C 279, the spectral slopes are 1.47 and 1.13 
respectively. That indicates the optical spectra of 3C 454.3 
is likely dominated by the synchrotron component. This in- 
terpretation is supported by the facts that 3C 454.3 strong 
variability in our observations and was in very high state 
(Gaur et al. 2012a). But 3C 279 seems to be in intermediate 
state during our observations and thermal component is also 
playing role with synchrotron component. 

5.4.2 ISPs 

For the BL Lacs showing intermediate behaviour, wc expect 
the aave Ssl. We have 11 ISP sources, and the average of 
their spectral indices is 1.19. This supports the hypothesis 
that synchrotron emission dominates the optical band of 
these sources too. 

5.4.3 HSPs 

The spectral slopes of the HSPs scatters in the range be- 
tween 0.5 and 1.8, while the predicted spectral slope under 
simple pure synchrotron emission models should be less than 
or equal to 0.75 or 0.80, corresponding to the spectral index 
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of optically thin synchrotron emitting plasma, athin (Urry 
& Padovani 1995). This indicates that the optical emission 
of HSPs is contaminated by other components, such as a 
thermal contribution from the accretion disc or host galaxy, 
or else includes significant amounts of non-thermal emis- 
sion originating from different regions of the relativistic jets. 
Also, HSPs are more generally stable in the optical bands, 
while great flaxes and SED changes axe frequently seen to 
be confined to the X- and 7-rays. This is consistent with the 
optical emission coming from a different region of the jet or 
with it being strongly diluted by thermal emission from the 
galaxy. For the nearby blazars, the host galaxy contribution 
is probably important (Pian ct al. 1994). For Mrk 421, Tosti 
et al. (1998a) determined the spectral slope to be a=0.85 
after the subtraction of an estimated host galaxy contribu- 
tion. We determined the spectral indices of 6 HSPs and for 3 
of them a could be found with the host galaxy contribution 
included as well as after the host galaxy subtraction (Table 
7). For Mrk 421, a is found to be 1.08 and 0.619 before and 
after the host galaxy contribution, respectively and for lES 
1959-1-650, a is found to be 1.03 and 0.67 before and after 
the host galaxy contribution is removed, respectively. Here 
we can say that there certainly is large deformation of the 
optical spectra of these HSPs due to the thermal contribu- 
tion from the host galaxies. However, for lES 2344+514, a 
isn't affected much, with high values of 1.77 and 1.61 be- 
fore and after the galaxy subtraction. So for this source, the 
optical emission may only be contaminated by non-thermal 
emission originating from elsewhere in the relativistic jet. 



5.4-4 Spectral slope variability 

We have investigated the possible correlations between the 
source intensity and the slope of the optical spectral en- 
ergy distribution. Fig. 7 shows the variability relations be- 
tween the spectral slopes and the R band fluxes (in mjy) 
for 12 blaaars. A linear regression analysis was applied to 
each blazar and the linear regression slopes b, their standard 
deviations, corresponding correlation coefficients and their 
respective p-values are listed in Table 7. Out of 15 sources 
studied, lES 0806-^524, lES 1218-^304, lES 1426-1-428, lES 
1959+650, 3C 454.3 and lES 2344+514 showed rather sig- 
nificant correlations (r > 0.6) between the source intensity 
and spectral slopes. From the above sources, a tendency of 
spectral flattening when brightening has been seen for HSPs. 
This trend has been seen for both HSPs and LSPs by other 
investigators (e.g., Spada et al. 2001; Sikora ct al. 2001; Chi- 
ang & Bottcher 2002; Vagnetti et al. 2003). However, the 
FSRQ 3C 454.3 has a tendency to show a spectrum that be- 
comes steeper when the object becomes brighter, which has 
also been seen by previous authors (Miller 1981; Ramirez et 
al. 2004). On the other hand, the FSRQ 3C 279 has shown 
the opposite tendency but it is quite weak, having r — 0.42. 
For the other 3 LSPs (3C 66A, S5 0716+714 and OJ 287), 
we found no significant correlation between flux and spec- 
tral slope. For J 0211+1051, ON 231, and BL Lac, we found 
negative correlations between spectral slope and R band flux 
(r >0.40) but they are not significant at a 0.95 confidence 
level. Only Mrk 421 showed a positive correlation between 
a and R magnitudes, but it is very weak, having r = 0.4. 



6 CONCLUSIONS 

We have carried out multi-band optical photometry of a 
sample of 11 blazars including 4 HSPs, 4 LSPs and 3 
ISPs during 2009 to 2011 on short-term timescales. We 
found significant flux variations in all the sources except 
for lES 1218+304 and lES 1553+113 on short-term time 
scales. Also, we searched for colour variations in the sources; 
only three sources (AO 0235+164, S5 0716+714 and lES 
1426+428) showed significant colour variations on short- 
term timescales. We have studied IDV for 4 sources (S5 
0716+714, OJ 287, lES 1426+428 and lES 1553+113) dur- 
ing 18 nights and found genuine IDV in 6 nights. We ex- 
amined blazar colours for IDV and found significant (B— R) 
colour variations on only 2 nights. 

We scached for possible correlations between colour and 
magnitude and found that 6 out of 10 BL Lacs followed the 
bluer-when-brighter trend; however, 5 sources (4 BL Lacs 
and 1 FSRQ) lacked any such correlations. We studied the 
optical spectra of a total of 15 sources (including 6 HSPs, 6 
LSPs and 3 ISPs) and compared these subclasses of blazars 
to provide some useful information on the emission compo- 
nents and emission mechanisms. LSP objects are extremely 
variable in optical bands and the average of the spectral 
slope of LSPs Q ~1.5, is in agreement with the synchrotron 
sclf-Compton model. However, the SED of FSRQs proba- 
bly have a contribution from a thermal bump that can no- 
ticeably contaminate the synchrotron component. The HSP 
spectral index varied from 0.5 to 1.8, which indicate that 
in addition to the synchrotron emission they have contribu- 
tions by components such as quasi-thermal emission from 
the accretion disc or non-thermal emission originating from 
different regions in the jet. For some HSPs, the thermal con- 
tribution of the host-galaxy plays an important role. This is 
clearly seen in Table 7 where we have been able to correct 
the spectral index for the host galaxy contribution for some 
nearby HSPs for which the galaxy could be resolved. The 
optical spectra of some of the sources become flatter when 
they brighten; however, for some sources we found no cor- 
relations or very weak opposite correlations. The relation 
between spectral slope and R magnitude of FSRQs suggests 
that the emission of FSRQs probably contains a thermal 
bump contribution and that this thermal component must 
be incorporated in detailed studies of individual sources. 
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Table 1. Details of telescope and instrument 



Site: 


NAO Rozhcn 


Telescope: 


60-cm Casscgrain 


CCD model: 


FLI P109000 


Chip size: 


3056 X 3056 pixels 


Pixel size: 


12 X 12 ;um 


Scale": 


0.330" /pixel 


Field: 


16.8' X 16.8' 


Gain: 


1.0 e-/ADU 


Read Out Noise: 


8.5 rms 


Binning used: 


2x2 


Typical seeing: 


1.5"to 3.5" 



^ With a binning factor of 1 X 1 
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Table 2. Observation log of optical photometric observations 



Source 


Q;2000.0 


(52000.0 


Date of Observation 


Telescope 


Data Points 


(z) 


(hh mm ss) 


(hh mm ss) 


(vyvy mm dd) 




(B, V, R, I) 


3C 66A 


02 22 39.61 


+43 02 07.80 


2009 11 13 


D 


2,2,2,2 


(0.444) 






2009 11 14 


D 


2,2,2,2 








2009 11 18 


C 


1,2,2,2 








2009 11 19 


C 


2,2,2,2 








2009 11 20 


C 


2,2,2,3 








2009 11 21 


C 


2,2,2,2 








2009 11 22 


C 


5,2,2,2 








2009 11 23 


C 


4,2,2,4 








2009 11 24 





3,3,2,2 








2009 11 25 


B 


2,3,2,2 








2009 12 21 


B 


2,2,3,2 








2010 11 01 


C 


2,2,2,2 








2010 11 03 


D 


0,4,4,4 








2010 11 04 


D 


0,2,2,2 








2010 11 05 


D 


3,2,2,2 








2010 11 06 


D 


4,4,4,4 


J0211+1051 


02 11 03.20 


+10 51 32.00 


2011 01 25 


D 


0,3,3,3 


(0.20) 






2011 01 26 


D 


3,3,3,3 








2011 01 27 


D 


2,2,2,2 








2011 02 07 


C 


2,2,2,2 


AO 0235+164 


02 38 39.93 


+16 36 59.27 


2009 11 13 


D 


2,2,2,2 


(0.94) 






2009 11 14 


D 


2,2,2,2 








2009 11 18 


C 


0,0,1,1 








2009 11 19 


C 


0,0,2,3 








2009 11 20 


C 


0,0,2,2 








2009 11 21 


C 


0,0,2,2 








2009 11 25 


B 


2,3,2,3 








2010 11 06 


D 


0,0,1,1 


S5 0716+714 


07 21 53.45 


+71 20 36.35 


2009 11 14 


D 


2,2,2,2 


(0.31+0.08) 






2009 11 15 


C 


2,2,2,2 








2009 11 16 


C 


2,2,2,2 








2009 11 19 


C 


2,2,2,2 








2009 11 21 


C 


2,2,2,2 








2009 11 22 


C 


2,2,2,2 








2009 11 26 


B 


2,2,2,2 








2009 12 11 


A 


0,0,61,0 








2009 12 13 


A 


1,1,105,1 








2009 12 15 


A 


1,1,115,1 








2009 12 21 


B 


3 j3^3;3 








2010 01 11 


C 


2,2,2,2 








2010 02 03 


C 


2,2,2,2 








2010 11 02 


C 


2,2,2,2 








2011 01 02 


C 


2,2,2,2 








2011 02 03 


C 


3,2,2,3 








2011 03 04 


c 


2,2,2,2 








2011 03 24 


c 


2,3,3,2 








2011 03 25 


c 


2,2,2,2 




ne no AQ 1 Q 






J) 


9 9 9 *> 


(0.138) 






2011 02 06 


c 


2,2,2,2 








2011 02 07 


c 


2,2,2,2 








2011 05 27 


D 


0,2,2,2 


OJ 287 


08 54 48.87 


+20 06 30.64 


2009 11 13 


D 


2,2,2,2 


(0.306) 






2009 11 15 


C 


2,2,2,2 








2009 11 20 


C 


2,2,2,2 








2009 11 21 


C 


2,2,2,2 








2009 11 21 


A 


1,1,1,1 








2009 11 22 


A 


1,1,35,1 








2009 11 22 


C 


2,2,2,2 








2009 11 26 


B 


2,2,2,2 








2009 12 21 


A 


1,1,132,1 








2010 01 10 


A 


1,1,90,1 
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Table 2. continued ... 



Mrk 421 11 04 27.20 +38 12 32.00 

(0.031) 



lES 1218+304 12 21 21.94 +30 10 37.11 

(0.182) 



ON 231 12 21 31.69 +28 13 58.50 

(0.102) 



3C 279 12 56 11.17 -05 47 21.52 

(0.5362) 



IBS 1426+428 14 28 36.60 +42 40 21.00 
(0.129) 



2010 


01 


11 


A 


1,1,80,1 


2010 


01 


11 


C 


2,2,2,2 


2010 


01 


20 


A 


1,1,95,1 


2010 


02 


03 


C 


3,3,3,3 


2010 


03 


16 


E 


3,3,3,3 


2010 


11 


02 


C 


2,2,2,2 


2009 


11 


15 


C 


2,2,2,2 


2009 


11 


21 


C 


2,2,2,2 


2009 


11 


22 


C 


2,2,2,2 


2009 


11 


26 


B 


2,2,2,2 


2010 


01 


11 


C 


2,2,2,2 


2010 


02 


03 


C 


2,2,2,2 


2010 


03 


06 


D 


2,2,2,2 


2010 


03 


16 


E 


2,2,2,2 


2010 


04 


03 


C 


2,2,2,2 


2010 


05 


14 


C 


2,2,2,2 


2010 


06 


07 


c 


2,2,2,2 


2010 


06 


10 


c 


2,2,2,2 


2010 


06 


12 


c 


2,2,2,2 


2011 


01 


02 


c 


2,2,2,2 


2011 


02 


05 


D 


2,2,2,2 


2011 


02 


08 


D 


2,2,2,2 


2011 


05 


02 


D 


2,2,2,2 


2011 


05 


02 


D 


0,2,2,2 


2011 


05 


27 


D 


0,2,2,2 


2011 


05 


29 


D 


0,2,2,2 


2011 


06 


23 


D 


0,2,2,2 


2011 


06 


26 


D 


0,2,2,2 


2011 


07 


06 


C 


2,2,2,2 


2009 


11 


15 


C 


2,2,2,2 


2009 


11 


21 


c 


2,2,2,2 


2009 


11 


22 


c 


2,2,2,2 


2009 


11 


26 


B 


2,2,2,2 


2010 


01 


11 


Bfr 


2,3,2,4 


2010 


02 


03 


C 


2,2,2,2 


2010 


03 


06 


D 


2,2,2,2 


2010 


03 


16 


E 


2,2,2,2 


2010 


04 


04 


C 


2,2,2,2 


2010 


05 


13 


C 


2,2,2,2 


2010 


05 


14 


C 


2,2,2,2 


2010 


06 


07 


C 


2,2,2,2 


2010 


06 


08 


c 


2,2,2,2 


2010 


06 


10 


c 


2,2,2,2 


2010 


06 


12 


c 


2,2,2,2 


2010 


08 


06 


c 


2,2,2,2 


2011 


02 


06 


c 


2,2,2,2 


2011 


05 


27 


D 


2,2,2,2 


2011 


05 


28 


D 


2,2,2,2 


2011 


05 


29 


D 


2,2,2,2 


2011 


05 


31 


D 


0,2,2,2 


2011 


06 


08 


B 


2,2,2,2 


2011 


06 


09 


C 


0,2,2,2 


2011 


06 


21 


c 


0,2,2,2 


2011 


06 


22 


c 


0,2,2,2 


2011 


06 


23 


c 


0,2,2,2 


2011 


06 


23 


D 


0,2,2,2 


2011 


06 


24 


c 


0,2,2,2 


2010 


03 


16 


E 


3,3,3,3 


2010 


03 


17 


E 


2,2,83,1 


2010 


03 


17 


A 


0,0,32,0 


2010 


04 


04 


C 


0,2,2,2 


2010 


05 


14 


C 


1,1,1,1 


2010 


06 


04 


c 


2,2,2,3 


2010 


06 


07 


c 


1,2,2,3 


2010 


06 


12 


c 


2,2,2,2 
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Table 2. continued ... 



lES 1553+113 15 55 43.04 +11 11 24.37 
(0.360) 



BL Lac 22 02 43.29 +42 16 39.98 

(0.069) 



2010 


06 


18 


F 


62,2,62,2 


2010 


06 


21 


F 


62,2,62,2 


2010 


07 


15 


D 


0,2,2,2 


2010 


07 


16 


D 


0,2,2,2 


2010 


08 


06 


C 


0,0,2,2 


2011 


02 


08 


C 


3,2,3,2 


2010 


03 


18 


E 


3,3,3,3 


2010 


05 


07 


C 


2,2,111,2 


2010 


05 


14 


C 


4,4,4,4 


2010 


06 


08 


c 


2,2,2,2 


2010 


06 


10 


c 


2,2,2,2 


2010 


06 


11 


c 


2,2,120,2 


2010 


06 


12 


D 


50,2,50,2 


2010 


06 


14 


D 


24,2,21,2 


2010 


06 


20 


F 


61,2,61,2 


2010 


06 


22 


F 


57,2,57,2 


2010 


07 


15 


D 


3,2,2,2 


2010 


07 


16 


D 


3,2,2,2 


2010 


07 


17 


D 


3,2,2,2 


2010 


07 


17 


C 


2,2,2,2 


2010 


08 


06 


C 


2,2,2,2 


2009 


07 


15 


C 


0,0,15,0 


2010 


06 


12 


c 


2,2,2,2 


2010 


06 


13 


c 


2,2,2,2 


2010 


06 


14 


D 


3,2,2,2 


2010 


06 


15 


D 


4,3,3,3 


2010 


06 


17 


D 


5,5,5,3 



A : 1.04 meter Sampuranand Telescope, ARIES, Nainital, India 

B ; 2-m Ritchcy-Chrcticii telescope at National Astronomical observatory Rozhen, Bulgaria 
Bfr : 2-m Ritchey-Chretien telescope with focal reducer (FoReRo2) at the National Astronomical Observatory Rozhen, Bulgaria 
C : 50/70-cm Schmidt telescope at National Astronomical Observatory, Rozhen, Bulgaria 

D : 60-cm Cassegrain telescope at Astronomical Observatory Belogradchik, Bulgaria 
E : 60-cm Cassegrain telescope at National Astronomical Observatory Rozhen, Bulgaria 
F ; 1.3-m Ritchey-Chretien telescope at Skinakas Observatory, University of Crete, Greece 
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Table 3. Results of intra-day variability observations 



Source Name 



Date 


Band 


N 


C-Tcst 

Cl,C2 


F 

Fi,F2,Fc(0.99),Fc(0.999) 


Variable 


A(%) 


11.12.09 


R 


64 


4.41, 


4.62 


19.44, 


21.36, 1.81, 2.21 


V 


12.89 


13.12.09 


R 


104 


2.07, 


1.75 


4.30, 


3.06, 1.59, 1.85 


PV 


7.71 


15.12.09 


R 


115 


4.85, 


4.61 


23.48, 


21.27, 1.55, 1.79 


V 


18.19 


22.11.09 


R 


35 


1.37, 


1.42 


1.87, 


2.01, 2.26, 2.98 


NV 


- 


21.12.09 


R 


129 


0.96, 


1.17 


0.92, 


1.36, 1.51, 1.73 


NV 


- 


10.01.10 


R 


90 


0.85, 


0.90 


0.73, 


0.81, 1.64, 1.94 


NV 


- 


11.01.10 


R 


80 


1.03, 


0.99 


1.06, 


0.98, 1.70, 2.02 


NV 


- 


20.01.10 


R 


62 


1.04, 


1.09 


1.09, 


1.20, 1.83, 2.24 


NV 


- 


15.03.10 


R 


31 


0.64, 


0.84 


0.40, 


0.71, 2.39, 3.22 


NV 


- 


17.03.10 


R 


83 


4.62, 


4.48 


21.42, 


20.09, 1.68, 2.00 


V 


11.99 


18.06.10 


B 


62 


2.62, 


2.34 


6.88, 


5.00, 1.83, 2.24 


V 


7.40 




R 


62 


2.00, 


1.99 


4.01, 


3.94, 1.83, 2.24 


PV 


3.69 




(B-R) 


62 


2.51, 


2.23 


6.29, 


4.96, 1.83, 2.24 


V 


7.70 


21.06.10 


B 


47 


3.06, 


2.80 


9.38, 


7.82, 2.01, 2.54 


V 


8.10 




R 


59 


3.14, 


2.66 


9.86, 


7.09, 1.86, 2.28 


V 


4.10 




(B-R) 


47 


3.33, 


3.17 


11.11, 


10.03, 2.01, 2.54 


V 


9.90 


07.05.10 


R 


111 


0.51, 


0.91 


0.26, 


0.83, 1.56, 1.81 


NV 


- 


11.06.10 


R 


120 


0.78, 


0.67 


0.60, 


0.45, 1.54, 1.74 


NV 




12.06.10 


B 


50 


0.57, 


1.00 


0.33, 


1.00, 1.96, 2.46 


NV 






R 


50 


0.47, 


0.80 


0.22, 


0.64, 1.96, 2.46 


NV 






(B-R) 


50 


0.54, 


0.92 


0.30, 


0.84, 1.96, 2.46 


NV 




14.06.10 


B 


24 


0.97, 


0.48 


0.94, 


0.23, 2.72, 3.85 


NV 






R 


21 


1.01, 


0.84 


1.02, 


0.70, 2.94, 4.29 


NV 






(B-R) 


21 


0.96, 


0.49 


0.92, 


0.24, 2.94, 4.29 


NV 




20.06.10 


B 


61 


0.70, 


0.78 


0.49, 


0.61, 1.84, 2.25 


NV 






R 


61 


0.78, 


0.90 


0.62, 


0.82, 1.84, 2.25 


NV 






(B-R) 


61 


0.67, 


0.86 


0.44, 


0.74, 1.84, 2.25 


NV 




21.06.10 


B 


57 


0.88, 


1.02 


0.78, 


1.04, 1.88, 2.32 


NV 






R 


57 


0.70, 


0.70 


0.49, 


0.49, 1.88, 2.32 


NV 






(B-R) 


57 


0.81, 


0.93 


0.66, 


0.86, 1.88, 2.32 


NV 





S5 0716-f714 



OJ 287 



lES 1426-f428 



V: Variable; NV: Non Variable; PV: Possible Variable 



Flux and spectral variations in blazars 17 



Table 4. Results of short-term variability observations 



Source Name 


Band 


N 


0-Test 
Ci, C2 


F 

Fi,F2,Fc(0.99),Fc(0.999) 


Variable 


A(%) 




B 


13 


2 80 2 Q5 


7 81 8 71 4 16 7 no 

I.O-L, 0.(i, rt.-LVJj 1 . \JyJ 


Y 


53.90 




Y 


15 


'i. -LO 5 0. yo 


1 7 91 1 5 R9 "^70 5 Q'^ 
-Li.z±, ±o.uz, o./u, o. yo 


Y 


43.51 




R 


15 


O.UO, I.UU 


95 4fi 91 1Q "^ 70 5 0*^ 
ZrO.iu, z_L.±y, o./u, o. yo 


Y 


43.50 




J 


15 




'i9 'iS 97 75 "^ 70 5 0*^ 
oz.oo, z/.(o, o./u, o. yo 


Y 


46.98 




yo- V ; 


12 


7.25 7.73 


59 R*^ 5Q f^9 4 4R 7 7fi 


Y 


26.10 




{D-l) 


12 


1 7S 9 '^S 


'^17 5 fi7 /] /Ifi 7 7f, 

O. -L 1 , O.U 1 , '^i.'l.yj. 1 . / U 


NV 






v-rtj 


15 


1.68, 1.77 


9 8*^ 14 S 70 5 0*^ 
z .oo, o.x^, o./u, o.yo 


NV 








1 ^ 

J-O 


1 99 70 


1 4Q fi^ ^ 70 5 


NV 










49 41 4"^ 1 Q 


1 7Q Q4 1 8fi 1 1 Q7 9Q 75 


Y 


92.52 




V 
V 


a 
u 


( .00, ( .00 


fi9 1 57 1 1 10 Q7 9Q 75 
uz.io, oi.i-L, ±u.y(, zy. ( o 


V 


4Q fio 
'iy .uu 




R 


14 


3.68 3.51 


1 5fi 1 2 "^l Q1 6 41 


Y 


57.30 




J 


16 


4.21. 4.32 


1 7 75 1 fi5 '^ 59 5 54 

-L ( . 1 O, J-O.UO, O.OZ, O.OI 


Y 


79.23 




("R Vl 




Q '^S 1 1 Q 


07 1 0*^ f^l 10 Q7 9Q 75 


Y 


79.54 








fi 1 n 97 


^^7 1 fi 97 75 1 Q7 9Q 75 
oi.xu, zi.io, ±u.y(, zy.io 


Y 


71.25 










Q fi'i fi 9'i 10 Q7 9Q 75 
y.uo, u.zo, J-U.y/, zy.io 


Y 


58.18 






14 




1 80 "iO 1 S4 7"^ Q1 fi 41 


V 


70 Q'i 
/ U.yo 




B 


18 


14 "^9 14 fin 


910 8^ 91^10 94 4 Q9 
z J.U.OO) z ±o. xy, O.Z1, i.yz 


Y 


122.96 




V 
V 


1 s 

J-O 


1 7 QR 1 7 Q7 
J- 1 .yu, J. 1 .y 1 


^99 fiQ "^99 QQ 94 4 Q9 
ozz.uy, ozz.yy, o.z4, 4.yz 


V 
V 


1 1 fi '^4 




R 


18 


9 9^ 1 HS 


5 04 S 59 S 94 4 Q9 

iJ.Ut:, O.O^, tj.iirtj '-t.cr^ 


Y 


112.87 




I 


18 


1 3 59 1 3 50 


1 84 63 1 82 1 4 3 24 4 92 


Y 


110.59 






18 


2 92 3 05 


8 54 9 30 3 24 4 92 

O.O^, ti/.OV./, O.^TT* 


Y 


14.25 






18 


R 1 S '^ 
U.XO5 u.uu 


"^S 91 '^1 SQ 94 4 Q9 


Y 


31.96 






18 


9 Qf; 1 Qfi 


5 52 1 85 3 24 4 92 


Y 


75.79 






12 


5.23 4.72 


97 '^S 99 '^9 4 4fi 7 7fi 

ZI.OO, ZZ.OZ, 'rt.'lU, i.lU 


Y 


21.38 


OJ 287 


B 


22 


QQ (S 40 


35 88 40 98 2 86 4 13 


Y 


117.31 




Y 


22 


^9 ^7 1 9 os 

J.Z. J. 1 , J.Z.tJO 


14S17 145 Q5 9 Sfi 41^ 


Y 


107.94 




R 


22 


4 89 4 8Q 


9*^ 99 9S Q9 9 8fi 4 1 S 


Y 


105.77 




I 


21 


R 41 R 5Q 


70 73 73 74 2 94 4 29 


Y 


102.89 




("R-V"; 

(DM) 


20 


0.55 1.47 


30 2 17 3 03 4 47 

V.OO, ii . i 1 , O.OO, '±.'±1 


NV 








20 


1.09 1.74 


1 ^^Q '^09 O'^ 4 47 
± .oy , o . uz , o. uo, '-t.'-ti 


NV 






CV-Vi) 


22 


1.11 0.55 


1 94 n '^n 2 8fi 4 1 S 

J_.ijrr, U.OO, ^.OO, rr.iO 


NV 








21 


U . \jO 5 U . OO 


4*^ ^^7 9 Q4 4 9Q 
u.io, u.oi, z.yi, 'i.zy 


NV 




Mrlc 491 


B 


29 


1.75 1.35 


^ 05 1 ^^9 9 4fi ^ 'ifi 
O .uo , J- .OZ , Z.'iU, o . ou 


NV 






Y 


29 


S 1 S OR 


fi5 fiK fi4 Q9 9 4fi '\ "ifi 
uo.uo, ut.yz, z.^u, o.ou 


Y 


95.38 




rv 


oo 


7 91 R Q4 


51 4S 15 9 ^9 ^ OQ 


V 


Q1 98 
y -L.zo 




J 


35 


4 9"^ 4 "=10 


1>^ 04 90 9Q 9 9fi 9 QS 


Y 


81.60 






9Q 


1 09 17 


1 04 0*^ 9 4fi ^ 
i.U'ri, U.UO, z.-tIu, o.ou 


NV 






CR ^^ 


29 


1.04 0.37 


1 OQ 1 "i 9 4fi 'i ^fi 
± .uy , u. xo, z.4:Uj o.ou 


NV 






(V R~l 


9Q 


1 19 7Q 

J. . J.Z, u. / y 


1 94 fil 9 4fi 'X 
± .z^, u.u 1 , z.^u, o.ou 


NV 






rR-T"! 

ysx i; 


33 


Q7 54 


n Q5 n 'HO 2 "^2 OQ 

u.yo> u.^^u, Zi.ijzi, Kj.Kjij 


NV 




1 1 91 S ^^04 


Y 




1 9R RS 

J-.ZrU, U.UO 


1 58 4fi 10 Q7 9Q 75 
1 . oo, u.'iu, J-U. y ( , zy . ( o 


NV 






R 




0.75 1.16 


5fi 1 '^4 10 Q7 9Q 75 
v.'.ou, -L.oi, J-U.y(, zy.fo 


NV 










U.UO, ±.\J^ 


1 00 OQ 10 Q7 9Q 75 
±.uu, u.uy, lu.yi, zy.io 


NV 








g 


n Qs "^0 

u.yo, u.ou 


1 00 OQ 10 Q7 9Q 75 
± . uu , u.uy, iu.y/, zy.io 


NV 








g 


1.78 1.17 


\ X7 10 Q7 9Q 75 
0.J.O5 x.oij ±u.yi, zy./o 


NV 




ON 231 


B 


34 


9 57 9 45 


fi fi9 5 9 9Q S 04 


Y 


65.64 




Y 


33 




8 00 8 97 9 "^9 "i OQ 
o.uu, o.zi , z.oz, o.uy 


Y 


133.47 




R 


35 


9 R4 9 70 


fi Qfi 7 9Q 9 9fi 9 Q8 
u.yu, ( .zy, z.zu, z.yo 


Y 


173.44 




J 


37 


7.74 7.80 


5Q 8Q fiO S4 9 90 9 KQ 
oy.oy, uu.oi, z.zu, z.oy 


Y 


159.48 




CR V) 


32 


QO R5 
u.yu, U.UO 


^^9 4*^ 9 '^5 "^15 

U.OZ, U.10, Z.OO, O.-LO 


NV 






(B-I) 


32 


1.19, 1.10 


1.41, 1.21, 2.35, 3.15 


NV 






(V-R) 


31 


0.79, 0.73 


0.63, 0.54, 2.39, 3.22 


NV 






(R-I) 


17 


1.60, 1.52 


2.57, 2.30, 3.37, 5.20 


NV 




30 279 


V 


11 


6.30, 5.81 


39.68, 33.75, 4.85, 8.75 


V 


43.84 




R 


11 


7.41, 7.66 


54.96, 58.73, 4.85, 8.75 


V 


55.86 




I 


11 


7.80, 7.38 


60.83, 54.45, 4.85, 8.75 


V 


41.81 




(V-R) 


11 


0.94, 0.99 


0.88, 0.98, 4.85, 8.75 


NV 






(R-I) 


11 


1.28, 1.41 


1.64, 1.98, 4.85, 8.75 


NV 




lES 1426+428 


B 


15 


5.26, 5.48 


27.70, 30.05, 3.70, 5.93 


V 


60.89 




V 


24 


6.68, 7.09 


44.67, 50.23, 2.72, 3.85 


V 


57.96 




R 


27 


5.73, 5.73 


32.88, 32.84, 2.55, 3.53 


V 


36.99 




I 


31 


4.08, 4.46 


16.63, 19.87, 2.39, 3.22 


V 


28.97 




(B-V) 


15 


2.63, 3.13 


6.93, 9.81, 3.70, 5.93 


V 


45.83 
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Table 4. continued ... 





(B-I) 


15 


4.17, 4.37 


17.36, 19.14, 3.70, 5.93 


V 


50.87 




(V-R) 


24 


3.06, 3.27 


9.38, 10.68, 2.72, 3.85 


V 


23.88 




(R-I) 


24 


3.89, 3.70 


15.17, 13.66, 2.72, 3.85 


V 


25.94 


lES 1553+113 


B 


36 


0.65, 1.05 


0.43, 1.11, 2.23, 2.93 


NV 


~ 




V 


32 


1.74, 0.75 


3.01, 0.57, 2.35, 3.15 


NV 






R 


34 


1.38, 1.15 


1.91, 1.32, 2.29, 3.04 


NV 






I 


30 


1.98, 1.96 


3.94, 3.83, 2.42, 3.29 


NV 






(B-R) 


31 


0.48, 1.08 


0.23, 1.18, 2.39, 3.22 


NV 






(V-I) 


28 


1.73, 0.73 


2.98, 0.54, 2.51, 3.44 


NV 




BL Lac 


B 


16 


3.46, 3.41 


11.96, 11.61, 3.52, 5.54 


V 


51.33 




V 


14 


16.90, 16.68 


285.44, 278.34, 3.91, 6.41 


V 


50.46 




R 


16 


42.04, 41.95 


1767.69, 1760.03, 3.52, 5.54 


V 


119.04 




I 


12 


25.92, 25.78 


672.10, 664.80, 4.46, 7.76 


V 


43.51 




(B-V) 


14 


0.76, 0.76 


0.57, 0.58, 3.91, 6.41 


NV 






(B-I) 


12 


0.57, 0.70 


0.32, 0.49, 4.46, 7.76 


NV 






(V-R) 


14 


1.40, 0.90 


1.96, 0.81, 3.91, 6.41 


NV 






(R-I) 


12 


1.79, 1.73 


3.20, 2.99, 4.46, 7.76 


NV 





V: Variable; NV: Non Variable 



Table 5. Fits to colour-magnitude dependences and colour- magnitude correlation coefficients on intra-day time-scale. 



Source name 


Date 




(B- 
m" 


-R) Vs R 

pa 


lES 1426-1-428 


18.06.10 


-0, 


,755 


13.473 






--0, 


,342 


0.0065 




21.06.10 


-0, 


,821 


14.529 






-0, 


,422 


0.0039 


lES 1553-1-113 


12.06.10 


-0. 


,425 


6.642 






-0. 


,180 


0.211 




14.06.10 


-1. 


,510 


21.069 






-0. 


,281 


0.231 




20.06.10 


-0, 


,402 


6.238 






-0, 


,228 


0.078 




22.06.10 


-0, 


,898 


12.812 






-0, 


,622 


3.894C-06 



m = slope and c = intercept of CI against V; r = Pearson coefficient, p = null hypothesis probability 
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Table 6. Fits to colour-magnitude dependences and colour-magnitude correlation coefficients on short time scale. 



Source name 


(B-V) 


vs V 


(V-R) 


vs V 


(R- 


-I) vs V 




(B-I) 


vs V 








m 


c 


m 


c 




m 


c 




j.a 


pa 


r 


p 


r 


p 




r 


p 


OKj \J\jr\. 


—0.131 


4.689 


0.012 


2.199 


—0.068 


2.856 


_Q 


.141 


3.412 




—0.233 


0.467 


0.030 


0.917 


—0.450 


0.092 




.199 


0.535 




n 1 


Q 709 




-LU.OUO 


—0 222 


/ .ooo 


n 

U 


A 70 


— u.ooo 




—0.081 


0.879 


0.778 


0.068 


—0.218 


0.678 


Q 


.322 


0.534 




0.055 


2.517 


0.006 


2.643 


—0.083 


3.104 


Q 


.046 


0.744 




0.639 


0.004 


0.070 


0.804 


—0.341 


0.213 


Q 


.183 


0.468 


OJ 287 


0.081 


2.190 


0.018 


2.401 


0.024 


1.758 





.106 


—0.004 




0.545 


0.016 


0.207 


0.381 


0.205 


0.386 





.418 


0.075 


Mrk 421 


0.137 


-0.541 


0.031 


-0.100 


0.317 


-3.065 





.137 


-0.541 




0.454 


0.013 


0.337 


0.0736 


0.872 


6.925C-10 





.454 


0.013 


lES 1218+304 


* 


* 


-0.089 


1.872 


0.167 


-1.832 




* 


* 








-0.400 


0.223 


0.590 


0.057 






* 


ON 231 


0.030 


1.921 


0.066 


1.010 


0.054 


0.911 





.079 


0.225 




0.001 


0.461 


0.440 


0.013 


0.518 


0.002 


0, 


.374 


0.029 


3C 279 


* 


* 


0.071 


-0.358 


0.203 


-2.770 




* 


* 




* 


* 


0.193 


0.714 


0.44 


0.383 




* 


* 


lES 1426+428 


0.467 


-6.694 


0.2357 


-3.38 


0.176 


-2.37 





.821 


-11.502 




0.386 


0.156 


0.577 


0.0031 


0.386 


0.076 





.569 


0.027 


BL Lac 


-0.021 


3.322 


0.050 


2.090 


0.026 


2.220 


0, 


.069 


1.434 




-0.131 


0.655 


0.436 


0.119 


0.345 


0.272 


0, 


.461 


0.131 



" m = slope and c = intercept of CI against V; r = Pearson coefficient, p = null hypothesis probability 

* missing entry is due to la^k of data 



Table 7. Flux density and spectral slope in the optical region. 



Source 


Class 


B [m 


■Jy] 


V [m 


■Jy] 


R [mjy] 


I [m 


Jy] 


a 




Ma 


b(+a) 


r 


P 






avg 


N 


avg 


N 


avg 


N 


avg 


N 


avg 


N 










3C 66A 


ISP 


8.13 


13 


11.42 


15 


13.02 


15 


15.41 


15 


1.01+0.15 


15 


0.48 


0.007+0.022 


0.087 


0.756 


J0211+1051 


LSP 


8.57 


3 


9.89 


4 


11.24 


4 


15.78 


4 


1.13+0.12 


4 


0.27 


-0.051+0.049 


-0.585 


0.415 


AO 0235+164 


LSP 


0.13 


2 


0.17 


2 


0.24 


2 


0.42 


2 


2.08+0.62 


2 


* 




* 


* 


S5 0716+714 


ISP 


9.73 


18 


13.06 


18 


15.92 


18 


19.34 


18 


1.20+0.16 


18 


0.52 


-0.003+0.006 


-0.112 


0.657 


lES 0806+524 


HSP 


4.05 


3 


4.61 


4 


5.10 


4 


5.79 


4 


0.63+0.02 


4 


0.05 


-0.130+0.032 


-0.945 


0.054 


OJ 287 


LSP 


4.47 


10 


6.59 


10 


7.56 


10 


10.16 


10 


1.41+0.14 


10 


0.44 


-0.016+0.024 


-0.229 


0.524 


Mrk 421 


HSP 


15.38 


29 


21.31 


29 


25.15 


29 


30.86 


29 


1.08+0.12 


29 


0.39 


0.010+0.007 


0.417 


0.156 






14.05 


29 


16.63 


29 


18.97 


29 


* 


* 


0.62+0.15 


29 


0.52 


0.011+0.007 


0.416 


0.156 


lES 1218+304 


HSP 


* 


* 


1.16 


6 


1.19 


6 


2.31 


6 


0.89+0.14 


6 


0.38 


-0.356+0.216 


-0.635 


0.175 


ON 231 


ISP 


2.57 


16 


3.32 


16 


4.26 


16 


5.37 


16 


1.30+0.14 


16 


0.53 


-0.083+0.041 


-0.474 


0.064 


3C 279 


LSP*" 


2.26 


4 


3.53 


11 


4.70 


11 


6.88 


11 


1.13+0.14 


11 


0.48 


-0.080+0.060 


-0.418 


0.200 


lES 1426+428 


HSP 




* 


1.053 


10 


1.38 


10 


1.74 


10 


1.36+0.21 


10 


0.73 


-0.893+0.390 


-0.624 


0.051 








* 


0.572 


3 


0.641 


3 


0.650 


3 


0.33+0.14 


3 


* 


* 


* 
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Figure 1. STV plots of blazars 3C 66A, AO 0235+164, S5 0716+714, OJ 287, Mrk 421 and lES 1218+304. For each source, the lower 
four panels show the calibrated LCs in B, V, R and I bands (upper LC), plotted with the differential instrumental magnitudes of standard 
stars with arbitary offsets (lower LC). The corresponding colour LCs are plotted in upper four panels in (B-V), (V-R), (R-I) and (B-t). 
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Figure 2. IDV light curves (upper panel is calibrated LC of blazar and lower panel is differential magnitude of standard stars) of S5 
0716+714 and lES 1426+428 along with their SF and DCF (from left to right) using R band data in each of the LCs. 
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Figure 3. IDV light curves of the blazars OJ 287, lES 1553+113 and lES 1426+426. X axis is JD(2455000+) and Y axis is Magnitude 
in each figure. In each panel, upper curve (triangles) is the calibrated light curve of the blazar and the lower curve (open circles) is the 
differential light curve of the standard stars. 
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Figure 4. As in Fig. 1 for STV plots of ON 231, 3C 279, lES 1426+428, lES 1553+11.3, and BL Lac. The lower right panel shows the 
colour variation (B— R) versus R magnitude on intra-day timescales with the source and respective date noted in each subpanel). 
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Figure 5. Color magnitude plots of blazars on short-term timescales. X axis is colour and Y axis is V magnitude in each figure. For 
lES 1553+113, Y axis is R magnitude. 
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Figure 6. Spectral indices, a against log(i/peafe) for LSPs, ISPs and HSPs. Red symbols are a's calculated from our data and black 
symbols are o's taken from literature. Triangles represent LSPs, squares ISPs and filled circles, HSPs. Open red circles represent the a's 
of four sources (lES 2344+514, Mrk 421, lES 1959+650 and lES 1426+428, from left to right) after host galaxy subtraction. 
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Figure 7. Dependence of the spectral slope on intensity for blaaars. X axis is spectral slope and Y axis is flux in R band (m Jy) . 



